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We present analytic computations of gauge invariant quantities for a point mass in a circular 
orbit around a Schwarzschild black hole, giving results up to 15.5 post-Newtonian order in this 
paper and up to 21.5 post-Newtonian order in an online repository. Our calculation is based on 
the functional series method of Mano, Suzuki and Takasugi (MST) and a recent series of results 
by Bini and Damour. We develop an optimised method for generating post-Newtonian expansions 
of the MST series, enabling significantly faster computations. We also clarify the structure of the 
expansions for large values of t, and in doing so develop an efficient new method for generating the 
MST renormalised angular momentum, v. 


I. INTRODUCTION 

The equations of motion of an extreme mass ratio bi¬ 
nary system — for example, a solar mass black hole or 
neutron star in orbit around a massive black hole 
are well approximated by a perturbative expansion of 
Einstein’s equations. With q = m/M (the ratio of the 
smaller mass to the the larger one) as an expansion pa¬ 
rameter, at zero-th order in q one obtains geodesic mo¬ 
tion in the spacetime of the larger mass. The order q 
corrections are commonly referred to as the first order 
self-force and are obtained by solving the linearized Ein¬ 
stein equations around the background spacetime of the 
larger mass, M. 

An alternative approach to the two-body problem - 
valid when the constituents are far apart — is the post- 
Newtonian approximation, which expands the Einstein 
equations in v 2 /c 2 , where v is a representative velocity 
and c is the speed of light. In the context of binary 
systems, the post-Newtonian expansion maps onto an 
expansion in 1/r, where r is the separation of the two 
objects. 

Recent years have seen synergistic development of the 
post-Newtonian and perturbative approaches to the ex¬ 
treme mass ratio binary problem. Comparisons of gauge 
invariant quantities computed in the two approaches have 
gone beyond simple (but valuable) crosschecks [lj. They 
have enabled impressive developments which would have 
been difficult, if not impossible, to achieve in each theory 
on its own. For example, on the post-Newtonian end self¬ 
force results have produced predictions of previously un¬ 
known terms in the post-Newtonian expansion [2|-[8(. On 
the self-force end, comparisons with post-Newtonian cal¬ 
culations and with Numerical Relativity simulations have 
given insight into the higher-order perturbative correc¬ 
tions not included in the self-force approximation @,E3. 

Traditional approaches to the perturbative treatment 
of Einstein’s equations have relied on numerical solutions 
of differential equations as a means to obtaining the lin¬ 
earized metric perturbation. More recently, alternative, 
functional methods have emerged as a compelling ap¬ 


proach to obtaining the linearized metric perturbation. 
In particular, there has been a surge of interest in ap¬ 
proaches pioneered by the work of Mano, Suzuki and 
Takasugi (MST) El , in which one writes solutions of 
the linearized Einstein equations in terms of a rapidly 
convergent series of hypergeometric functions. 

Functional methods have the distinct advantage of giv¬ 
ing a representation in terms of exact quantities, rather 
than a truncated numerical value. Since modern com¬ 
puter algebra systems can efficiently evaluate hyperge¬ 
ometric functions to an essentially arbitrary number of 
digits, a numerical approach built on top of functional 
methods is a powerful tool. Indeed, a parallel pair 
of works by Shah and Pound 11211 . and by Johnson- 
McDaniel, Shah and Whiting |13( has used numerical 
functional methods to simultaneously obtain many of the 
same results as we present here. Reassuringly, a direct 
comparison of our results with those of [12|,[l3| has shown 
the two to be in perfect agreement [T3 | , providing a strong 
independent check to both calculations. 

In this work, we compute post-Newtonian expansions 
of the linearized metric pertubation for a circular-orbit 
binary system, ignoring spin terms, up to order y 21 , 
where 1/y = (Mff) _2//3 is an invariant measure of the 
radius defined through the orbital frequency, fl. From 
this, we compute very high order post-Newtonian ap¬ 
proximations of all known gauge invariant quantities up 
to quadrupole order. Our approach works in the Regge- 
Wheeler gauge (with a transformation to the asymptot¬ 
ically flat post-Newtonian gauge) and is fundamentally 
based on functional methods, but avoids any numerical 
evaluation (and the associated numerical truncation). In¬ 
stead, our method builds on a series of developments by 
Bini and Damour EHm Much of the calculation de¬ 
scribed here relies heavily on the methods they devel¬ 
oped, combined with some modifications which allow the 
calculation to be efficiently taken to much higher post- 
Newtonian order. 

The layout of the paper is as follows. In Sec. [TT] we 
give details of our method, including an MST-based ex¬ 
pansion for low Umultipole modes and an exact large-£ 
approach for the higher multipole modes. In Sec. IIIII we 
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give explicit expressions for the metric perturbation in 
Regge-Wheeler gauge in terms of homogeneous solutions 
of the Regge-Wheeler equation. In Sec. IIVI we summarise 
our results and we conclude with a discussion in Sec. m 
Finally, in Appendix [A] we give equations for the gauge 
invariant quantities we compute (written in terms of this 
metric perturbation) and in Appendices iBl IQ iDl lEl FI 
and Owe give our high-order post-Newtonian expansion 
of the gauge invariant quantities. As the results become 
increasingly unwieldy with increasing post-Newtonian or¬ 
der, we restrict our printed results to order y 15 - 5 and opt 
instead to provide the higher-order terms electronically 

III- 

This paper follows the conventions of Misner, Thorne 
and Wheeler [20j; a “mostly positive” metric signature, 
(—,+,+,+), is used for the spacetime metric, the con¬ 
nection coefficients are defined by T x v = \g XrT {g a ^, v + 
9crv,p - g^,a), the Riemann tensor is R a \^ = T^ - 
+ ~ r^r^, the Ricci tensor and scalar are 

R a p = R^aiifi and R = R a a , and the Einstein equa¬ 
tions are G a p = R a p — ^g a pR = 87 t T a p. Standard 
geometrized units are used, with c = G = 1, but we 
include the explicit dependence on G and c in our post- 
Newtonian expansions in cases where they are convenient 
for post-Netwonian order counting. We use the spherical 
coordinates {i, r, 9, <f>} for the background Schwarzschild 
spacetime and write tensors in terms of these coordinate 
components. 


II. METHODS 


with 


F(a, b , c, C) = r( :?? ) 2 Fi(a, b, c; C), 


F(c) 


x = 1 — 


2 GMuj 


2GM’ 


However, we find it more useful to instead work with the 
following decomposition: 


2 F 1 (a,b, c; C) = 
r(c)r(6 — a) 


r(6)r(c — a) 


(l-C) a 2 F 1 (a,c-b,a-b+l, J ^) 


r(C)r( °—S(! - 0 ~ b 2 F 1 (c - a, b, b - a + 1 , 


T(a)r(c — b) 


( 2 . 2 ) 


Thus F (a, b, c, £) = F\ (a, b , c, Cf) + F 2 (a, b 1 c, £) with 
F 1 (a, 6, c, C) = 

r(a)T(b- a) _ c) -a 2 j p i( c _ b a _ b + ! 

1 (c — a) ** 

(2.3) 

F 2 (a,b, c, C) = 

T(b)r(a - b ) 


F(c — b) 


(l-C) b 2 Fi{c — a, 6, 6 — a + 1, yz^)- 


(2.4) 


In our calculation of the homogeneous solutions to 
the Regge-Wheeler equation, we assume much the same 
strategy as [161] in approaching the various l mode 
regimes. Namely we use the analytical results of Zerilli 
| 21 j with asymptotically flat gauge correction for £ = 0 , 1 , 
generate the analytical MST series solutions for ‘low’ t 
modes, and use a variation of Bini & Damour’s post- 
Newtonian ansatz valid for large-1' values to tackle the 
rest. 


A. Low i solutions 

The MST expansion of X^(r) 

The homogeneous solutions of the Regge-Wheeler 
equation are given in [ll| as series of hypergeometric 2 F 1 
functions and irregular confluent hypergeometric U func¬ 
tions. Following a similar notation to [16:], the horizon 
solution may be written as 

OO 

XZ(r) = Cfa(x) J2 a > 


As a note, this leads to the X^ n = Xq + X 0 " 1 repre¬ 
sentation of Eq (2.15) in lllf if we let n —> —n in the F\ 
sum. Once again following |T(| we will keep track of the 
large-r, small-w coupled expansions by using the small 
parameter 77 = 1/c, and the Bini & Damour variables 
X\ = GM/r , A 2 1 / 2 = wr. In this way, each instance of 
Xi and X 2 must each come with an rj 2 . 

At this point we can do an examination of the lead¬ 
ing powers of 77 against n for the various terms in the 
sum. The structure of the MST series coefficients a n 
is described in Sec. 6 of (llj, but we present it again 
in terms of powers of 77 and in a form (Table [T| that 
makes it straightforward to combine with other elements 
of the expansions. We note that the corresponding be¬ 
haviour is also discussed in the Living Review by Sasaki 
and Tagoshi [22] for the MST series coefficients for the 
Teukolsky function, which simply differ from those of the 
spin-2 Regge-Wheeler equation by a complex conjuga¬ 
tion fp . 

This is not the only irregular behaviour we need to 
account for; the dependence on 77 of both the parameters 
and argument of the 2 -F 1 also go through order jumps, 
which are tabulated in Table [II] In constructing this and 
subsequent Tables the key elements are that 


F(n + v — 1 — ie, —n — — 2 — ze, 1 — 2 ie, x), 

( 2 . 1 ) 


Au = v — i • 


1 — x 
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\n < -21- l|-2l < n < -£ - l\n = -i\ n = -£+ 1 \n > -£ + 2 
a n | 3(|w| — 1) | 3(|n| + 1) | 31 + 6 13(1 - 1) + 61 3|n| 


TABLE I. The leading behaviour of the MST coefficients for the spin-2 Regge-Wheeler equation in terms of powers of r/. For 
example, a„ = 0{rf^) for n > — l + 2. 


One subtlety here as that while we tend to think of 2 -F 1 
as being regular in its argument, this is only true if the 
parameters are held fixed. It is easy to see, for exam¬ 
ple, that the expansion of the hypergeometric function 
in Eq. (12.31) is given by 

! (n-\-£— 1 — ie-\-Ais)(n-\-£+3—ie+Av) v ^2 | 
i_1 2(n-\-£+l—Av) “• ’ 

and the denominator is manifestly rf when n + £ +1 = 0 , 
leading to a divergent limit as 77 —> 0. This divergence 
is more than compensated for by the behaviour of the 
prefactor, but leads to larger term than a casual analysis 
would suggest. This term corresponds to the daggered 
entries in Table m 

Combining the behaviour of the various elements gives 
a master table, Table m for the orders of each term in 
the series. This turns out to be the most useful table 
for optimising our calculations of the inner solution since 
they allow us to say with certainty what we do and do 
not need to calculate to correctly give a desired order. 
For example, one can see immediately from looking at 
these that for n > 0 our contribution from F i will die 
out quickly and can be largely ignored. 


The MST expansion of X^(r) 


The solution satisfying the boundary condition at in¬ 
finity is 


+oo 


X£(r)=Cf up) (z) E <(~2iz) n x 

n=—oo 

U(n + v + 1 — ie, 2n + 2v + 2, — 2iz), (2.5) 

with 


C{ up) (z) =e iz z v+1 { 1 - l)-^ 2 !/ e - 7 re e- l7r(l,+1) 


( 2 . 6 ) 

(2.7) 


and z = In a similar vein to the inner case we split 
this solution by making use of the identity 

r (6 — l) 


+ 


T(a) 


z °M (a — b + 1, 2 — 6 , z). 

( 2 . 8 ) 


Thus U (a, b, () = U\ (a, b, () + C / 2 (a, b, () 


Ui (a, 6,C) 


U 2 (a, 6,C) 


T(a — 2)T(a)T (1 — 6) 

= r(a* + 2)T(a*)r(a-6+l) 


M (a, 6, C) 


(2.9) 


r( n — 2 )r (6 — i) i~b 

T(a* + 2 )T(a*) 

M (a- 6+ 1,2 -6,0- (2-10) 


Analysis of the behaviour of the M functions 
for small 77 produces Table IIVI Again we must 
take due care that parameters as well as the ar¬ 
gument depends on 77 . For example, when consid¬ 
ering Ui(n + £ + l + Av — ie, 2(n + £ + 1) + 2Ais, —2 iz) 
the Taylor coefficients in the expansion of M are regular if 
n+£ > 0 but if n+£ = — k ( k £ N) the Taylor coefficients 
of order (— 2iz) 2k ~ 1 and beyond are order 7 y~ 3 . This 
means that for k = 1 the function behaves as rj~ 2 , while 
for k > 2 the behaviour is regular but we need to include 
3 more terms than one might have expected to get to 
the appropriate 77 order. Correspondingly when consid¬ 
ering U 2 (n + £ + 1 + A is — ie, 2(n + £ + 1) + 2Av, —2 iz) 
the Taylor coefficients in the expansion of M are regular 
if n + £ < 0 but if n + £ = k (t £ N U {0}) the Taylor co¬ 
efficients of order (—2 iz) 2k+1 and beyond are order r]~ 3 . 
This means that for k = 0 the function behaves as 77 ” 2 , 
while for k > 1 the behaviour is regular but we need to 
include 3 more terms than one might have expected to 
get to the appropriate 77 order. 

Combining the behaviour of the various elements gives 
a master table, Table El for the orders of each term in 
the series for X up . 


B. Large £ solutions 

As described in [l(|, one can generate ‘large-£’ homo¬ 
geneous solutions using the template 

X“( PN) =r t+1 [1 + 7 fA™’ e + 77 *A: N ’ e 

+ r, e A™’ e + ...] 

x up( p N) =r _ f + + ^ A PN,-e-i 

+ r, e Aj; N ’- e - 1 + ...] (2.11) 

where the Ai are for the most part polynomials in X\ and 
A 2 1 / 2 , but one finds extra r dependent log terms appear¬ 
ing after the sixth order, log 2 terms appearing after the 
twelfth order, etc. However, the precise details appear 
opaque. 
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n < -£- 3 

n = -l - 2 

n = —£ — 1 

n = —£ 

n = — l + 1 

n > -l + 2 

Ei 

2n + 21 — 5 

-9 

— 19 T 

-11 

-9 

2n + 2t - 5 

f 2 

-2 n -2£-7 

-9 

-11 

- 11 + 

-9 

—2 n -21-1 


TABLE II. The leading behaviour of the hypergeometric functions appearing in Eqs. (f2~3l) and (12.41 . The terms marked with 
a t highlight the interplay between parameters and argument; see the text for a full discussion. 



n < —2i — 1 

-21 < n < -£ - 3 

n = -e- 2 

n = — £ — 1 

n = — l 

+ 

1 

II 

s 

n> -£ +2 

y M+7 a„Fi 

r? l+7 a n F 2 

\n\ + 4Z — 1 
5|n| — 3 

|n| + 4£ + 5 

5|n| + 3 

6£ + 7 
6£+7 

5£ — 6 4 

5£ + 2 

5^ + 2 
5£ + 2+ 

5^+1 

5^+1 

3|n| + 2n + 4l + 2 
3|n| — 2 n 


TABLE III. The combined behaviour following from Tables HI and HT1 


In generating the MST expansions for high values of £ 
and for general t above the given order the full structure 
of this expansion becomes clear: 

(1) In the PN expansion of the MST series using the 
Bini-Damour variables Ai, X 2 the combination 

e = 2 XiX 2 1/ V 

is r-independent. A phase of the form 

oo oo 

M /up = ^ yui/up e i = J- ^ /up (2X 1 X 2 1 / 2 r j 3 y 
2 — 0 2=0 


greatly simplified the algebra. However, for large-f we 
may take an alternative approach. Acting on a function 
/(Ai, A 2 1(/2 ), we have 


df _ 

dr yVl dXi 


+ A , 1 / 2 


df 

dA 2 V 2 ’ 


and use this substitution in the spin-s Regge-Wheeler 
operator, which is given by 


(1 - 2X 1 7f) 2 r-^~ 
dr 



+ (l-2X lV 2 )(4X lV 2 



+ r) 2 X 2 - (1 - 2X lV 2 ) (£{£ +!) + (! — s 2 )2Xnf) 


is therefore irrelevant as it will drop out of the Green 
function via normalisation through the Wronskian. 

(2) Having factored out an r-independent phase of the 
type described in (1), the MST expansions for X' n may 
be expressed in the form 


X 


in(MST) _ £ 1 5 l °(«,2n( 2 - YlA ' 2 1/ V) 2j 

tm — e -'M X 

[l + ? 7 2 A 2 + rf A\ + rf Aq + ...] 


( 2 . 12 ) 

where the Ai are now strictly polynomials in X\. X 2 
which we may express as 


A e 2 n =Y!^ 2nti) X i^A, 172 ) 2 ^-^ 
2=0 
n 

= J2' a {2n,i)XyX 2 n -\ 

i —0 


Requiring that the ansatz (12.121) satisfy the spin-2 Regge- 
Wheeler equation immediately determines all coefficients 
at orders upto n = 2£. This means that, as Bini & 
Damour did, we may determine at any given order the 
necessary expansion without solving the full MST expan¬ 
sions. (As a check we have solved the MST expansion for 
large £ explicitly to order 77 21 and found full agreement 
with the much simpler method described here.) 

An additional bonus of this method is that one may 
identify the coefficients d(6j, 2 j) as being precisely the co¬ 
efficients of the renormalised angular momentum v_ cor¬ 
responding to —l — 1 , that is 

OO 

Me) = -l- 1 (2T3) 

1 =i 

or, in other words, we have 

win(MST) v-{2X x X, 2 1/ V) v 

"Wm “ e A 1 X 

[l + V 2 Ai + ^A{ + V 6 Ai + ...]. (2.14) 


where the prime indicates that we omit terms where the 
combination is a power of X\X 2 ! 2 as these are already 
included in the phase term. This occurs if * = 2(n — i) 
which is only possible if 2 n = 6 j and i = 2j for some 
j £ N, explaining our notation for the phase expansion. 

In our calculations in the previous subsection we cal¬ 
culated Aj n from the hypergeometric series to order 77 40 
up to £ = 20 explicitly, and then rewrote it in the form 
( 12 . 121 ) prior to calculating the metric perturbation, as it 


(3) The full MST expansions for A up are in general much 
more complicated, involving, for example, high powers of 
logA 2 . Nonetheless, the leading terms are again easy; 
with the ansatz 


A, 


up(MST) _ x/2 1 £M’2T( 2A 'i a 2 V V) 2j 

lm 6 V^2 ) 


[l + r, 2 B'+M + - + /4 + 0 (m, 

(2.15) 
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n < -i- 3 

n = -£- 2 

n = —l — 1 

n = —t 

n = —l + 1 

n > — i + 2 

Ui 

0 

-3 

—5 T 

-6 

-6 

-3 

u 2 

-2 n -U-A 

-3 

-5 

—6 f 

-6 

—2n -2£-l 


TABLE IV. The leading behaviour of the hypergeometric functions appearing in Eqs. (EH) and (I2T01 .The terms marked with 
a t highlight the interplay between parameters and argument in this case. 



n < —2£ — 1 

—2t < n < —E — 3 

n = -l- 2 

n = — l — 1 

n = —t 

n = — l + 1 

n > -£ + 2 

V™ +L a n (-2iz) n Ui 

V 2t+1 a n (-2iz) n U 2 

2|n| +2^-2 

4|n| — 6 

2|n| +21 + 4 
4|n| 

At+ 5 

At+ 5 

At + F 

4^+1 

At+ 1 
41! + l f 

4t - 1 
At- 1 

3|n| + n + 2£ — 2 
3|n| — n 


TABLE V. The combined behaviour following from Tables HI and IIVI 


the Regge-Wheeler equation again determines all the 
given coefficients (that is, orders up to n = 21) and 


v+(e) — i + ^ &( 6 j, 2 j) e2j - (2-16) 

i=i 


In both cases it is easy to see why we can determine to 
this order and no higher since the r independent factor 
( 2 X 1 X 2 1 / 2 ry 3 ) ± i 2Z+1 ) will serve to mix them at that order. 
To obtain the correct solutions we must supplement the 
Regge-Wheeler equation with the boundary condition; 
our simple ansatz cannot implement these while the MST 
expansion does. 


III. RECONSTRUCTED METRIC 
PERTURBATION IN REGGE-WHEELER GAUGE 

In Regge-Wheeler gauge, the components of the met¬ 
ric perturbation may be written in terms of the homoge¬ 
neous solutions, X?' m and X??. For the case of a point 
mass on a circular orbit of radius ro, the metric per¬ 
turbation components are given explicitly below for the 
case r < ro (we ignore terms involving S(r — ro) and 
its derivatives since they are not relevant to our calcula¬ 
tions). Here, W = X^(r 0 )X^'(r 0 ) - X^(r 0 )X^(r 0 ) 
is the Wronskian, A = (£ — l)(i'+2)/2 and A = A+ 3 M/r. 
To obtain the components for r > ro, simply interchange 
X?? and X?? in these expressions. 

The (t, t) component of the metric perturbation is 
given by 


rw ^(00, <h)Yim(8, m^XZ(r) + h%XZ'(r))(h^XZ(r 0 ) + K? X^/M) 
8A(A + l)r 3 ro /2 W(r 0 - 2¥)(r 0 - 3Af) 1 / 2 (9m 2 M 3 + A 2 (A + l) 2 r 3 ) 


(3.1) 


where 


h\ n t =(6M - 2(A + l)r)(4Am 2 Mr 3 + 4A(A + l)rg(2M - r)), 
h% =2r(2 M - r)(12m 2 M 2 r 2 + 4A(A + l)rg(3Af - r)), 

h Y = - 8(6M 3 (—A + m 2 - 1) + 3M 2 r 0 (2A 2 + 3A - m 2 + 1) - A(A + l)Mr 0 2 (A + m 2 + 4) + A(A + 1 ) 2 rg), 
hif =8r 0 (2M - r 0 )((A + 1)(3 M 2 + Ar 0 (r 0 - 2 M)) - 3m 2 M 2 ). (3.2) 


The (t, r) component of the metric perturbation is given by 


h 


RW 

tr 


where 


M ztt mY; m (Q 0 , ^(h^X^jr) + ^^'(r))(Cg(>^) + K?'X%'(r 0 )) 

r 0 - 3 M 4A(A + l)r 2 r4W(2M - r)(2M - r 0 )(9m 2 M 3 + A 2 (A + l) 2 r 3 ) 


(3.3) 


/i£ =12m 2 M 2 r 4 - 4rg(18M 3 - 3(2A + 5)M 2 r - 3(A 2 - 1 )Mr 2 + A(A + l)r 3 ), 
h£ = - 4A r 2 rl(r - 2Af)(-3Af + Ar + r), 

hi? =8(6M 3 (—A + m 2 - 1) + 3M 2 r 0 (2A 2 + 3A - m 2 + 1) - A(A + l)Mr 2 (A + m 2 + 4) + A(A + l) 2 rg), 
hi?' = - 8r 0 (2M - r 0 )((A + 1)(3 M 2 + Ar 0 (r 0 - 2 M)) - 3m 2 M 2 ). (3.4) 
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The (r, r) component of the metric perturbation is given by 

h RW = nY; m (9o,MY tm (0, + h^X^'(r))(h^XZ M + Kg'X?*' (To)) 

8A(A + l)rr 5 0 /2 W(r - 2M) 2 (2M - r 0 )\Ar 0 - 3M(9m 2 M 3 + A 2 (A + 1 ) 2 rg) 

where 


/C =(6M - 2(A + l)r)(4Aro 2 Mr 3 + 4A(A + l)rg(2 M - r)), 

=2r(2M - r)(12m 2 M 2 r 2 + 4A(A + l)rg(3Af - r)), 

= — 8(6M 3 (—A + m 2 — 1) + 3M 2 r 0 (2A 2 + 3A — m 2 + 1) — A(A + 1)Mto(A + m 2 + 4) + A(A + l) 2 rg), 

K p' =8r 0 (2M - r 0 )((A + 1)(3M 2 + Ar 0 (r 0 - 2 M)) - 3 m 2 M 2 ). (3.6) 


The (t, 0) and (t, (f>) components of the metric perturbation are given by 


where 


H, = 


Ke = -imHt esc 0Y £m (6,4>), h™ = H t sin 8Y lm f(6, </>), 

~KT~ 4tt Y; m j(e 0 , fo)(2M - r)(rl“'(r) + XZ(r))(r 0 XZ'(r 0 ) + ^M) 


ro — 3M A(A + l)rroVF 

The (r, 0) and (r, <f>) components of the metric perturbation are given by 

h™ = —imH r esc 0Y) m ($, <j>), h™ = H r s'm 6Y im j(0, 0), 


where 

H 4inY; m j(e 0 ^ 0 ) m Mr 2 XZ(r)(r 0 XZ\r 0 ) + (r 0 )) 

A(A + l)ro /2 W(2M - r)(r 0 - 3M) 1 / 2 
The ( 6 , 0) and (<(>, <(>) components of the metric perturbation are given by 

"■ee — A, — sin OK, 

where 

K ■xY? m {0 o ,<MYi m {0, <t>){K*X? m {r) + {r)){K**XZ M + 

8A(A + 1) r r 5 0 /2 W(r 0 - 2 M)(r 0 - 3My/ 2 (9m 2 M 3 + A 2 (A + l) 2 r 3 ) 

with 


(3.7) 

(3.8) 

(3.9) 

(3.10) 

(3.11) 

(3.12) 


K ia =24m 2 M 2 r 3 - 8(A + l)rg(-6M 2 + 3 Mr + A(A + l)r 2 ), 

K in ' = - 8(A + l)rr 3 (r - 2M){3M + Ar), 

X up =8(6M 3 (—A + m 2 - 1) + 3M 2 r 0 (2A 2 + 3A - m 2 + 1) - A(A + l)Mr 2 (A + m 2 + 4) + A(A + l) 2 r^), 

X up ' = - 8r 0 (2M - r 0 )((A + 1)(3M 2 + Ar 0 (r 0 - 2M)) - 3m 2 M 2 ). (3.13) 


The ( 0 , <(>) component is exactly zero, as required by the Regge-Wheeler gauge conditions. 

I- 


The above expressions may only be used for the mul¬ 
tipole modes £ > 2. For completeness, we also give the 
metric perturbation components for l = 0,1, which were 
derived analytically by Zerilli 211. As in Ref. (l6j , a cor¬ 
rection term is added to these to Zerilli’s solutions to 
account for the change to an asymptotically flat gauge 


compatible with the one used in post-Newtonian theory: 


=o 2^1-f 2E X 

2 M H ( r o -r) + ——H(r - r 0 ), 


// =u - 
n tt — 


r o 


hi=° = 


ro 1 - 

2 Ei 




>H{r - r 0 ), 


(3.14) 
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2 M 
r 0 


with El = ^rrm- 

from the odd sector, 


For £ = 1 we have a contribution 



ro) , 

(3.15) 


and a contribution from the even sector, 


1 - 


sin 9 cos <j>H(ro — r), 


/ 4 =1 = - 2E 
u r(r - 2 M) 

h^ 1 = 6EQ r f' 1a s ^ n ^ s ^ n ^^( r o “ r)> 




(r - 2M) 2 
r(ro — 2M) . 


^ = - 6E (r-2M)3 


sin0cos^Ff(ro — r), (3.16) 


where 17 is the Heaviside function and L\ = y m(t 0 - 3 M) • 


IV. RESULTS 

Our main result is the high-order post-Newtonian ex¬ 
pansion of gauge invariant quantities for a circular orbit 
extreme mass ratio binary system, ignoring contributions 
from the spin of the constituent bodies. Similar post- 
Newtonian expansions have been given to lower orders in 
other works [ 3 , H, [H, El E1EH Isll - IH ; we have verified 
that our results are in perfect agreement with these. 

Our post-Newtonian expressions for the invariants take 
a standard form in all cases, and is given as a series in 
powers of y and logy. For example, the spin-precession 
invariant is given by 

Ai/> = c 2 y 2 + c 3 y 3 + c 4 y 4 + (c 5 + cf log y) y 5 + 

+ (c 6 + Cg n log y) y 6 + c 6 . 5 y 6 ' 5 + (c 7 + c 1 " log y) y 7 
+ C 7.5 V 7 ' 5 + (c 8 + 4 n log y + cf log 2 y) y 8 ■ ■ ■ , 

(4.1) 

where we use the notation that c n corresponds to the 
coefficient of y n , and cjf is the coefficient of y n log (y) k . 
The coefficients are exact numerical quantities involving 
rational numebers, logs of integers, 7 r, Euler’s constant 
( 7 ), and the Riemann zeta function (() with integer ar¬ 
gument (these various numerical quantities can be seen to 
come about from the infinite sums over £ and m spherical- 
harmonic modes). We give the explicit values for many 
of the coefficients in the appendices, with higher-order 
coefficients being available electronically 0 - 

I 


In addition to gauge invariant quantities, our expres¬ 
sions may also be used to compute a high-order post- 
Newtonian expansion of the linear-in-mass-ratio piece of 
the effective one body radial interaction potential. As 
was shown in Refs. [29 h 31| . this has a straightforward 
expression in terms of the redshift invariant, 

a(y) = -(l-3y) 3 / 2 At/-^=^, (4.2) 

V 1 - 3y 

so substituting our post-Newtonian expansion for A U in 
to this expression yields and expansion for a(y). 


V. DISCUSSION 

Although we have given post-Newtonian expansions 
for all currently known invariants for circular orbits in 
Schwarzschild spacetime, there remain further quantites 
one could compute. For example, a forthcoming work 
[ 32 } will present numerical and post-Newtonian approx¬ 
imations for octupolar invariants for circular orbits in 
Schwarzschild spacetime. 

From an astrophysical perspective a particularly com¬ 
pelling direction for future study is the extension to the 
Kerr background corresponding to a rotating black hole. 
The derivation of the various gauge invariants for circular 
orbits in Kerr spacetime has already been worked out in 
[H and the MST method in general applies equally well 
to the Kerr case. The outstanding issue is therefore the 
development of efficient post-Newtonian expansion tech¬ 
niques for the Kerr case. Fortunately, we anticipate that 
many of the tricks employed in this work will carry over 
to the Kerr case; we will explore this in more detail in a 
future work. 
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Appendix A: Gauge invariants 


In this work we evaluate the redshift, spin-precession and tidal-tensor invariants by writing them in terms of the 
components of h a b in Schwarzschild coordinates. In particular, for our particular choice of off-worldline extension the 
invariants are given by 
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12 
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12 ht 
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In the above we evaluate all of the metric perturbation components using the expressions given in Sec. ED with 

(M) = (f,o). 


Appendix B: Post-Newtonian expansion coefficients for AC/ 


In this Appendix, we give the analytic post-Newtonian expansion coefficients for Detweilcr’s redshift invariant [24j . 
AC/, up to order y 155 . We have also computed higher order coefficients up to order y 21 ’ 5 , but they are too long to 
give here; instead, we make them available electronically fl9| . 
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Appendix C: Post-Newtonian expansion coefficients for A ip 


In this Appendix, we give the analytic post-Newtonian expansion coefficients for the spin-precession invariant 
A ip, up to order y 155 . We have also computed higher order coefficients up to order y 21 ' 5 , but they are too long to 
give here; instead, we make them available electronically {Ts]. 
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216432216 lu gV z y LU &\°) 448238748672000 

_ 2950513264271989525498007 , 18487328 >/ Q \ , 122783649257283464 ^,_ 1380030544 2 

11420261125842318750 ' 1225 ^ 1491353238375 ' 385875 ' 1 231525 

+ PIP l°g 2 ( 2 ) - PIP l°g 2 ( 3 ) + 47 ^gs wa log(2) - ^ § ^ 1 °g( 2 ) 

+ 111 98i96098000 569 log(3) - ^7^(3) - ™ log(2) log(3) - ™ 0 4 7 4 ° 3 3 2 125 log(5), 

30579194748876386 __ 690015272 ,,_ 14320872 1/91 _ 1108809 , /o\ 

1491353238375 385875 ' 42875 ^g^/ 343 LU S>\°h 

_ 345007636 
1157625 > 

106979456 a/ o\ , 2084022912123714512565994459849 , 190712970681326946869 _ 5723400896 2 

23625 i 50082608917173919086000000 n ' 12655197479925000 ' 7r 2480625 ' 71 

I 177395403334811 3 , 106979456 . 3 , 999808 5 _ 160490433344 , 2/ 0 \ _ 2538986352360389949581 1 /q\ 

' r 93648555000 1 “ r 70875 ' n ' 10125 71 17364375 ' g v > 88586382359475000 1 ^gv ^/ 1 

53470676864. _,„„/ol , 499725952 _3 ,„„/r,5 , 33855365577141 _ ,,„/ol , 29176220703125 _ ,„„/cl 

- 5788125 77Tl0g(2) + 165375 7 T l0g(2) + 31 4 91 4 600 o 77 log( 3 ) + 2623 688496 ^^gl 5 ). 

190712970681326946869 _ 5723400896 . , 53489728 3 _ 26735338432 1 no ./o', 

25310394959850000 2480625 ’ n _t ' 70875 5788125 lu g^ z ^> 


In 2 _ _ 1430850224 
c 12.5 2480625 71 > 
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C 13 = 


62650400736283736 /-to\ 79334586176 .,a/p5 , 1209920 a/k1 , 189559817152 a/pM.„./o5 3947937408 /■/ o 

147496474125 W 1091475 7s('j) + 21 < 3V i: V+ 467775 13475 S W i,J B W 

55918518949617438556735975663717584848969 _ 31820709089016321579091597087 _ 230884798395817567208 

394794945767967357415052250000 32222266766564102353125 ' 2214659558986875 ' 

--- " 268326297213415329697937454609374713 " 

73789210202161733959680000 n 


r ln — 
c 13 — 


Jn 2 

c 13 

c !n 3 

c 13 

C13.5 


l 195609221811488 3 , 268326297213415329097937454609374713 2 _ 48902305452872 2 

t 11345882625 7 73789210202161733959680000 2269176525 l n 

, 19750691096766759888050303278403159 4 , 3037644237389156207824805 6 __ 5990948040507941 8 

_t " 37751549977222749093888000 n ~r 53198770598313984 77 549755813888 71 

_ 61457225667104 1 n „ 3 / 9 , , 155944342152 i no .3/o', , 242362228261971549865288 l no . 2 /o', _ 4035505222260896 -. 1 n r 2 /o', 

138929175 1U B 5187875 1U B ' 112947637508330625 1U B 11345882625 7 iu 6 \ L ) 

i 467833026456 i„„2/om (o\ 308594519722437276357 i „„ 2 /o\ i 467833026456 ., l 2/o\ 

5187875 i0 8 l z ; 10 gW 367253406520000 i0 & W 4 5187875 7 10 B W 

, 467833026456 1 ~„ 2 /q'i , 299718641201171875 l„„2/r\ , 14180966230199741 i„„2/ 7 \ 

+ 5187875 io gl z J 10 g W + 794739097152 10 § + 307181160000 10 § 

_ 205476770552879570318960688943 i /o', i 2333171461995431294768 i /o\ _ 73182403279456 .,2 1/91 
46952445859850549143125 lu gv z 7 4~ 3227075357380875 / lu Bv z / 1620840375 ' 1 U B'> Z 7 

I 18295600819864 2 1 /o', , 11090055998662708077^46193027 1 /o\ _ 308594519722437276357 1 / 0 \ 

' 972504225 lu Bv z J 4~ 15198908564002042880000 LU S>\°) 183626703260000 7 iu 6 l°i 

, 467833026456.,2 i„„/o\ 38986085538 _2 i„„/o\ 409010222195260682949 i„„/oM_ fo\ , 935666052912 i„„/oM„„/o'i 

+ 5187875 7 l 0 S( 3 )- 1037575 77 lo g( 3 ) - 183626703260000 - lo g( 2 ) lo g( 3 ) + 5187875 7 lo g( 2 ) lo g( 3 ) 

_ 933945255836406040673828125 1 /r', , 299718641201171875 _ 1 /cl , 299718641201171875 1/91 1 /r\ 

148007479166482231296 “r 397369548576 7 iu BW4" 397369548576 lu Bl z J 1 U BW,I 

_ 2447069882207180299994701 1/71 , 14180966230199741 . 1 /y\ , 14180966230199741 ] no ./ 7 \ 

19614031164389376000 153590580000 7 1 LI BV‘J ' 153590580000 iu Bl z J iu BV J; 

_ 39667293088 A/0'1 _ 31837810312556420711589661087 _ 231202291990193925608 , 97804610905744 .,2 

1091475 'sW 64444533533128204706250 2214659558986875 I' 3781960875 ' 

_ 24451152726436 2 _ 2017752611130448 l n „ 2 /o', , 233916513228 1 2/01 , 1164811253436376604824 1 /o', 

2269176525 n 11345882625 1U B \ > ' 5187875 u ° v 0 > ' 3227075357380875 iu gv z l 

73182403279456.,i„„/ol 308594519722437276357 i„„/ol 1 467833026456., 1 _ to\ 1 467833026456 i„„/oM„„/o5 

1620840375 7 10 Bl z j 367253406520000 10 6 W 4 5187875 7 i0 gW 9 5187875 10 £A Z J io gW 

1 299718641201171875 i /c\ 1 14180966230199741 1 / 7 , 

' 794739097152 9 307181160000 iu BVU) 

_ 38168822370617659002 , 48902305452872 . _ 18295600819864 1 /o', 

2214659558986875 ‘ r 3781960875 ' 1620840375 1 U BV Z J 


. 4451152726436 

11345882625 


116958256614 

5187875 


1166092488 a/ o', _ 1187374613661695186615691845565290057 , 245883449484395004266099 _ 751070271916 . 

42875 2281975776845387078391936000000 71 1612272158942445000 l w 121550625 ■ 

_ 13163514596884821173 ^3 , 388697496 . 3 , 562644568 5 , 1381319157332 1 2/ 0 x _ 28829034 1 2/o' 

178962388605000 ' ' 42875 > n ' 496125 1 ' 121550625 1 u o \ > 2401 1 u o 


r ln 

c 13.5 

„ln 2 

c 13.5 

C14 


1316351459688482U373 3 , 388697496 . 3 , 562644568 5 , 1381319157332 1 2/^1 _ 28829034 1 

178962388605000 ' ' 42875 > n ' 496125 “ ' 121550625 U B \ > 2401 h L 

3081012717124731361878079 1/93 _ 214793576 . 1/93 , 48046504 3 1/93 , 2133354786803205 

8061360794712225000 n 1U BV Z J 324135 lu Bl z J ' 926 I 71 iu BV z ^ "r 982942940980C 

57658068. , 4435236 _3 i„„/ol 57658068 _ i„„/o\i /o\ 1688840218427734375 _ i„„/cl 

2401 77Tl0g(3) + —343—77 log(3)- 2m~ 77 l°g( 2 ) l°g( 3 ) - 11699027003664 ^ lo g( 5 ) 


SfM^logtS), 


1 ouoiuif 1 1 < iz-i <oijuio 1 ou 1 74 Intri 9 1 _ £± 2. 1 1 _ lntri 9 l I 

' 8061360794712225000 71 324135 AU &V Z / ' 9261 

- log(3) + log(3) - PIPtt log(2) log(3) _ 

_ 245883449484395004266099 _ 751070271916 . , 194348748 3 __ 107396788 l /o\ 

3224544317884890000 121550625 '' ' r 42875 324135 1U BV Z J 

_ 187767567979 
121550625 7r ’ 

_ 25105689326437633143222302545119678455060078765034473 , 5294485011363817767374 a/q\ _ 113555673327320608 a/o\ 
44918728131563797834340291702753325000000 1328795735392125 245827456875 

1998903808.,2 a/o1 , 499725952 _2 a/o\ , 18681344 a/oi2 89492704864 a/c\ 37362688 _,a/c1 349184a/t1 

- 55125 ^ C(3) + 33075 n C(3) + 5 25 C(3)- 70375 -C(5)-525-74(5)-~C(") 

76204544 a/oM _2/ol 142268027251758368 a/oM _/o\ 23998569472 .,a/oM„„/o\ , 238553295129 a/oM„„/o! 

- 525 -C(3) log (2)- 27314161875 -C(3) log(2)- 165375 74(3) log(2) + --C(3) log(3) 

1 943603515625 a/o\ _ 224285696 a/ r\ 1 / 9 \ , 28465651655704059830901242446584177479 . 

"T" 764478 ^OJlUg^o; 1575 l,yO) 4 - 4194696298784653172534930156250 7 

_ 149954059365283950864513030853 2 _ 14862394284976579472 .3 , 53470676864 .4 

144152246060944668421875 7 33219893384803125 7 17364375 7 

, 1382576599003842137026984607581609563901342613 2 , 3715598571244144868 2 _ 26735338432 .2 2 

' 79248535357017878238402183168000000 ' 6643978676960625 ' n 3472875 ' 1 

1 178409537584128378780724226998135283471761 4 , 499725952 4 , 1618696988208768975806817129793 6 

' 87066550718668735126230269952000 ” l " 826875 < n ' 7149914768413399449600 

_ 2547114068232098727607 8 , 856314857344 l n0 -4/o', , 151607779358849976638576 l r , 0 -3/o' ) 1 5137575533056 1 3/o', 

791648371998720 17364375 ^ 33219893384803125 ° 1 / ' 52093125 ' 0 \ > 

_ 33855365577141 i n0 .3/o\ _ 145881103515625 l S/ri _ 2277262395801345716119706739413 l no . 2 /o', 

787286500 io B W 1967766372 1U B 78254076433084248571875 iU g V z ) 

1267756872949042315472 1 2 / 0 x , 4076943104 .,2 1 2/ 0 x _ 1019235776 2 1 2 _ 101566096731423 1 2/r,N 1 /o', 

316379936998125 ’ 10 & V Z J 4 55 T 25 > i0 8 V Z J 33075 77 10 & V z ; 787286500 10 & l z J 10 gW 

75421247572218176728578441 l no .2/o\ _ 101566096731423 i„„2/o^ _ 101566096731423 1 /o\ l no . 2 /Q', 

33926135187504560000 iu g 787286500 ' lu & 787286500 1U BV Z ^ 1U B 

__ 145881103515625 l no . 2 /o') i n0 ./c'\ , 69917822132336320439453125 l no . 2 /c', _ 145881103515625 .,1 2/r\ 

655922124 iu g lu gv' J / ) 50877570963478267008 & 655922124 ' lu & 

__ 145881103515625 1 /o\ l n0 ' 2 /'c,'\ _ 4542009473315817564229 l no . 2 ('7' 1 1 86184302376552379542620466756226641423347 1 /o', 

655922124 iu Bv z J 1U B 5357853792720000 V ‘ 7 4" 1833082282568893436397764478281250 lu Bl z 7 

_ 41360868951644100025764392413246 . 1 /o\ , 10411999054067186375792 2 1 /o', , 1283923466752 3 1 /o', 

2738892675157948700015625 / 1U BV Z 7 “r 11073297794934375 < 1U B^ Z 7 ' 52093125 ' lu gv z 7 

_ 2602999763516796593948 _2 1 /o\ _ 320980866688 .,_2 1 /o\ , 2999821184 4 1 ../ 0 , 

6643978676960625 '* 1U BV Z 7 10418625 ' iu &V z 7 "f 2480625 H lu Bl z 7 

3837129087597300548202099012782283 1 /o\ 75421247572218176728578441 ., 1 /o', 101566096731423 .,2 1 /o', 

100288912523241479797760000 iu gl°^ 4- 16963067593752280000 ' 787286500 ' lu Bl°j 

4 . 33855365577141 2 1 /o', , 131731477811538857834428809 1 /o', , /o', _ 101566096731423 1 /o', loff fQ') 

+ 629829200 ~ 10 gv' 3 J + 16963067593752280000 10 Bl z J 10 BV ,3 J 393643250 7 io gl z ^ iOgl'Jj 


+ 
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r ln 

'-14 


r> ln 

-14 


r ln - 
c 14 — 

c 1 " 4 — 
c 14 — 

C14.5 = 


_ 14522555334465885013333427526900390625 i /r\ , 69917822132336320439453125 i /r\ _ 145881103515625 2 i /r\ 

1743425022700190479786142662656 1 U 6 W ' 25438785481739133504 / lu 6655922124 ' 1U E>( 0 J 

, 729405517578125 „2 i — , 69917822132336320439453125 i—/oM_145881103515625 .. i — com_ 

“i 7871065488 17 10 gW H 25438785481739133504 327961062 7 io gl z J 

, 2248625818037574093829749927087227 l / 7 \ _ 4542009473315817564229 l ( 7 \ 

' 131141373768223806873600000 lu gw 2678926896360000 

_ 4542009473315817564220 i l( 7 \ _ 81402749386839761113321 l /i i \ 

2678926896360000 lu gl‘J 2054356066344960000 iu gv 1 J 7; 

28462003489715167876845043737346977479 _ 56777836663660304 a/q\ _ 1998903808 a/ o\ _ 18681344 cn 
8389392597569306345069860312500 245827456875 ^7°' 55125 525 

_ 14999284736 a 10\1 /oi __ 149770950317944669377889830853 _ 7431197142488289736 .,2 , 106941353728 3 

165375 SW 1 U 6 '^/ 144152246060944668421875 ' 11073297794934375 ' ' 17364375 ' 

, 1857799285622072434 2 _ 26735338432 2 , 249862976 4 , 2568787766528 l n a- 3 ('9'| , 633878436474521157736 

' 6643978676960625 3472875 ' ' ' 826875 71 ‘ 1 “ 52093125 g V > 316379936998125 u ° v > 

, 4076943104 .. 1 2/q\ _ 101566096731423 1 2/o\ _ 145881103515625 l no .2/r\ _ 20678006470941745421284906926623 1 /o-i 

^ 55125 ' lu o 1574573000 iu g 1311844248 ° W 2738892675157948700015625 lu g^ z ' 

, 10411999054067186375792 i /q\ 1 641961733376 2 1/91 _ 160490433344 2 1 /q\ 

11073297794934375 / iu gl z J f 17364375 ' 10418625 iuf 7 z J 

, 75421247572218176728578441 1 /o\ _ 101566096731423 .. 1 lol _ 101566096731423 1 /o\ 1 lol 

33926135187504560000 iu gW 787286500 ' iu 6 >W 787286500 lu gW 

, 69917822132336320439453125 1 /r\ _ 145881103515625 .. 1 (r\ _ 145881103515625 i /n\ 1 /cl 

50877570963478267008 655922124 ' Lu &\°) 655922124 iu gW 

_ 4542009473315817564229 1 / 7 \ 

5357853792720000 lu 6 V ‘ A 

_ 150067369071353402594431110853 _ 499725952 f lo 1 _ 3715598571244144868 , 26735338432 2 _ 6683834608 2 

576608984243778673687500 55125 "s'. 0 / 11073297794934375 ' ' 5788125 ' 3472875 1 

, 1019235776 l„ n .2/o'i , 2602999763516796593948 , 320980866688 _ 101566096731423 

55125 1U g \ z > + 11073297794934375 ~r 17364375 7 10 gl z J 3149146000 

_ 145881103515625 1 /'r\ 

2623688496 iu gv J 7 

_ 1857799285622072434 , 26735338432 , 160490433344 1 /<j\ 

33219893384803125 ' 17364375 ' ' 52093125 lu 8\ z 7 

3341917304 

17364375 ’ 

_ 5309850881930378361213063000400507073797099 _ 31697702965724519600778001907 . 

9009966035282625018582291963648000000 115321796848755734737500 

_ 960372800793516616313 3 _ 871515123500656 . 3 , 125822671264 5 

71 8031832000592400 34037647875 ' ,l ' 147349125 


'364375 ’ 

871515123500656 aid! 
11345882625 71 *7 


r ln 

c 14.5 


r ln 

c 14.5 


iiO^OOOZDZO J - ' yUUyyOOlMOZOZOZDUlOOOZ 

2017887846520122448 2 _ 960372800793516616313 ^ _ », x,,xox^^ _ 

39313483295625 ' n 8031832000592400 n 34037647875 l /l ~~ r 147349125 ,l 

486929066270811952 1 nP -2/o\ , 166561812368496 i 2/o\ , 05792943231500199837619674037 i /q\ 
1572539331825 /l & ^ ^ 1997331875 & W ~ r 115321796848755734737500 iU &W 

_ 529898658421540832 i /q\ ■ 795346539665488 3 l nr /o^ _ 24616880040022202752230447 l nr /o\ 
7862696659125 ' t[ lu ol z / ' 6807529575 n 15623878046877100000 n 

+ r log(3) - 46 ™ 456 ^ 3 log(3) + 333 1 1 ^g ,a 7rlog(2) log(3) 

, 13713804^6935830078125 i n „/g\ , 42358546109606626367 
“ r 20286112824353376 71 7 509152772700000 lu gl'l; 

_ 31799589853782602615359096307 , 2017887846520122448 _ 435757561750328 3 

230643593697511469475000 39313483295625 34037647875 1 

_ 264949329210770416 i 166561812368496 i /o\ 

7862696659125 ‘ ' 1997331875 71 1U 6>V% 

504471961630030612 

3Q31 34839Q5R9K 71 ) 


Cl5 = 


_ 10732014117075454744748966 //o\ 

6387987609740269795047242089727919825000000 ' 2447198812680496875 

331375078473130816 .,^791 3353344 _,2 s/o\ i 4191680 i 1976320 , 550606278944 Alcl 3952640 

81942485625 9 7 SW i 27 71 S l' 3 - 1 7 3 606375 ‘’W 3 

2240000C(7) - M™92 c(3) log 2 (2) _ 159668496 C(3) log 2 (3) + 1132 g^S 1568 C(3) log(2) 

- 71 T2 8 8 3 625 576 7C(3) l0g(2) + ^ C( 3) log(3) - ®992 7C(3) log(3) _ 319336992 c(3) log(2) log(3) 

- 868 g 8 7 054 125 C(3) log(5) - » ((5) log(2) - a^iac(5) log(3) 

. 1444626847565977343847871701769271761693931. . 743974276294437883505577844651591 2 80654004101045801163952 3 

225502163588004442967953125 ' 143952871334146875 ' 

1675885675742282873^2 , 201635010252" 1 - 































































































































17 


345948408 i rin .2/ 9 \ i„_.2/o\ , 8444201092138671875 l„„2/ 9 \ i zr\ _ 18043699048009895547239775390625 1...2, 
2401 10 g 10 6 \°) ' 2924756750916 i0 g W lo &l 0 7 703738561566831389254656 10 ° 1 

8444201092138671875 ., l 2/r\ , 8444201092138671875 j z 9 \ l no -2zrc\ , 2124911309544887559242617 l nD .2z 7 \ 
2924756750916 ' iu & 2924756750916 1U B \*>) ~r 338787811021271040 ° V' / 

5914802257134088557178479343028167512162443 l Zo\ i 21582307719896314698020411808946354 l zrA 
288127206051056068321066798449843750 430504130486190300211546875 7 1U BW 

542228830585444590986576 .,2 i„ rr 7 9 \ , 12300964281344 


r ln — 

c 15 — 


i ZUDUOiUODUDM^iUDD /y5^yb<±L)i viu ' 

542228830585444590986576 2 i z 9 \ i 12300964281344 3 i z 9 \ , 135498760158289169941844 2 1 z«\ 

91606372667184375 ' g^ ' t 364651875 ' > 54963823600310625 g^ z 

3075241070336 . 2 i z 9 \ i 8988539072 4 i z 9 \ _ S212 132810613748046832911 8774308357987909 i Zo^ 

72930375 < n iu Bl z 7 ' 1929375 iu gV z / 888278500142161147403448064000000 lu gW 

37695796473493246706069582373039 i zo\ _ 6657887258847138663 2 i zo\ , 230632272 3 l /ol 

1493120051543155235200000 7 1 u b7°7 2233961229500 7 1 U 6 WJ + 2401 7 lu gW 

2219295752949046221 _2 i zo\ 288290340 . _2 i„„zo\ , 13305708 _4 „ to \ 

- ! 787168983600 - 77 lo g( 3 )-2401 - ^ log( 3 ) + m5 7T log( 3 ) 

518915434980801888780465029425029 i z 9 \ i zo\ _ 91708992333763253613 _ i z 9 \ i /ol 
16424320566974707587200000 lu Bl z j lu BV°j 12286786762250 7 lu gl z ) lu g7 0 7 

691896816 .,2 i_zoM_Zo\ 288290340_2 i„„zo\ 1_Zol i 335449641407501498004901489998569392578125 i„„ZE\ 

2401 ^ 1 ( Jgl z j 10 SW 2401 71 10 Bl“J iOgl'Sj “I 1507190932124314669775120331866112 10 gl°J 

18043699048009895547239775390625 l Zr\ , 8444201092138671875 2 i zr\ _ 42221005460693359375 2i zr\ 

351869280783415694627328 ' ' 2924756750916 ' u gv d 7 35097081010992 u gW 

18043699048009895547239775390625 i zo\ l no .Zc'| , 8444201092138671875 _,i z 9 \ l Zcl 

351869280783415694627328 iu gl z 7 1 U BW 1 " 1462378375458 ' lu gV z ^ lu gv J 7 

10864198815859511930364238491875158331 l Z 7 \ , 2124911309544887559242617 . i z 7 \ 

80979798301878200744448000000 lu Bl‘;~r 169393905510635520 ' g'' 1 ' 

21249113095448875592 12017 i z 9 \ l no .Z 7 \ , 379462616391725675509799117 l Zi i \ 

~ r 169393905510635520 iu Bl z ^ lu gV ') "r 333832860781056000000 iu gv 1J J! 

1443095427096991114707380066473385508493931 , 165174494992853408 /-{q\ _ 3353344 azq\ _ 1976320 z/r\ 

301847549196344452526831884090312500 ' 81942485625 9 'SW 3 

_ 35954156288 1 z 9 \ __ 159668496 AZol 1 zo\ , 744221366722721348748916146091591 _ 40327002050522900581976 . 

128625 SV°7 1 U BV Z / 343 SW iu gW~r 225502163588004442967953125 ' 47984290444715625 

1 644182528 3 , 10051135161735879628694 2 _ 161045632 2 , 419168 4 _ 3995392675328 l no -3z 9 \ , 115316136 l r , 0 -3zo\ 

' 19845 ' 28790574266829375 3969 '' 135 52093125 g 7 > 2401 g v 0 ' 

_ 9165637237632884112617272 1 2z 9 \ _ 580251275776 i„„2zoi , 345948408 l„„2z 9 \ 1 zo\ _ 6657887258847138663 l r , n .2zo'i 

1007670099339028125 1U g V z 7 13505625 7 lu g 7 Z 7 ~r 2401 lu g 4467922459000 1U g 
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Appendix D: Post-Newtonian expansion coefficients for AAi 


In this Appendix, we give the analytic post-Newtonian expansion coefficients for the tidal invariant AAi (see [8| for 
a definition) up to order y 15 - 5 . We have also computed higher order coefficients up to order y 21b , but they are too 
long to give here; instead, we make them available electronically 0- 
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Appendix E: Post-Newtonian expansion coefficients for AA2 


In this Appendix, we give the analytic post-Newtonian expansion coefficients for the tidal invariant AA 2 (see Q for 
a definition) up to order y 15 - 5 . We have also computed higher order coefficients up to order y 21 ' 5 , but they are too 
long to give here; instead, we make them available electronically 0- 
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Appendix F: Post-Newtonian expansion coefficients for AA 3 


In this Appendix, we give the analytic post-Newtonian expansion coefficients for the tidal invariant AA 3 (see @] for 
a definition) up to order y 15 - 5 . We have also computed higher order coefficients up to order y 21 - 5 , but they are too 
long to give here; instead, we make them available electronically 0- 
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Appendix G: Post-Newtonian expansion coefficients for AA fl 


In this Appendix, we give the analytic post-Newtonian expansion coefficients for the tidal invariant AA S (see @] 
for a definition) up to order y 15 - 5 . We have also computed higher order coefficients up to order y 21 ' 5 , but they are 
too long to give here; instead, we make them available electronically [li¬ 
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